Thermal Pressures in Neutral Clouds inside the Local Bubble, as 
Determined from C I Fine-Structure Excitations^ 

Edward B. Jenkins 

Princeton University Observatory 
Princeton, NJ 08544-1001 

ebj@astro.princeton.edu 
ABSTRACT 

High-resolution spectra covering the absorption features from interstellar C I 
were recorded for four early-type stars with spectrographs on the Hubble Space 
Telescope, in a program to measure the fine-structure excitation of this atom 
within neutral clouds inside or near the edge of the Local Bubble, a volume of 
hot (T ~ 10^ K) gas that emits soft x-rays and extends out to about 100 pc away 
from the Sun. The excited levels of C I are populated by collisions, and the ratio 
of excited atoms to those in the ground level give a measure of the local thermal 
pressure. Absorptions from the two lowest levels of C I were detected toward 
a Del and 6 Cyg, while only marginal indications of excited C I were obtained 
for 7 Ori, and A Lup. Along with temperature limits derived by other means, the 
C 1 fine-structure populations indicate that for the clouds in front of 7 Ori, 6 Cyg 
and a Del, 10^ < p/k < 10^cm~^K at about the ±lcr confidence level in each 
case. The results for A Lup are not as well constrained, but still consistent with 
the other three stars. The results indicate that the thermal pressures are below 
generally accepted estimates p/k > 10^ cm~^ K for the Local Bubble, based on the 
strength of x-ray and EUV emission from the hot gas. This inequality of pressure 
for these neutral clouds and their surroundings duplicates a condition that exists 
for the local, partly-ionized cloud that surrounds the Sun. An appendix in the 
paper describes a direct method for determining and eliminating small spectral 
artifacts arising from variations of detector sensitivity with position. 

Subject headings: ISM: atoms - ISM: lines and bands - techniques: spectroscopic 
- ultraviolet: ISM 



^ Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. 
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1. Introduction 

The Sun is immersed in two concentric volumes of interstellar material with very different 
properties. In the immediate surroundings [out to distances ranging from 0.05 to several pc 
away (Redfield & Linsky 2000)], the gas is partly ionized, with characteristic properties 
0.10 < n(H I) < 0.24 cm-3, 0.04 < n{e) < 0.10 cm-^ and T ^ 7500 K (Lallement 1998; 
Puyoo & Ben Jaffel 1998; Ferlet 1999; Gry & Jenkins 2001). This cloud, often called the 
Local Interstellar Cloud (LIC), is surrounded by a large cavity containing fully-ionized, very 
hot, low-density material with n(H+) ~ 0.008 cm~^ and T 10^ K (Berghofer et al. 1998; 
Burrows & Guo 1998) that is prominent in soft x-ray emission (Snowden et al. 1990, 1998). 
This large volume of hot gas is believed to be the product of probably several supernova 
explosions (Cox & Reynolds 1987; Breitschwerdt & Schmutzler 1994; Maiz-Apellaniz 2001; 
Smith & Cox 2001) and is known as the Local Bubble (LB). The LB is also conspicuous by 
showing a general absence of cold, neutral gas up to a well defined boundary (Sfeir et al. 
1999; Vergely et al. 2001), as is shown in Figure 1. 

A persistent puzzle has been the apparent mismatch in thermal pressures between the 
two media. Using the parameters stated above, one finds that 1400 < p/k < 3600cm"^K 
for the LIC, which contrasts with an apparent representative value of 16, 000 cm~^ K for the 
fully-ionized, hot gas in the LB. The objective of the research presented here is to help answer 
the question, "What thermal pressures are found for other neutral clouds within the Local 
Bubble? That is, do they have values similar to the LIC, or are they closer to matching 
the apparent pressure of the LB?" To gain insights on this question, one can observe stars 
located behind individual clouds inside the LB and measure their absorption features of C I 
arising from different fine-structure levels of the ground state. The ratios of populations 
of these states are governed by an equilibrium between coUisional excitations (governed by 
local densities and temperatures) and spontaneous radiative decays. 

Stars suitable for viewing the C I features had to satisfy a number of conditions to yield 
useful results. Their selection is described in §2; ultimately four such stars were observed 
in a manner described in §3. The C I absorption features are very weak, and in order to 
obtain useful measures of their strengths particular care was exercised to remove spurious 
signals arising from the detector, as outlined in §4.1. (Mathematical details about the 
correction method are presented separately in the Appendix.) Section 5 describes how the 
equivalent widths of various absorption features were combined and corrected for saturation 
(very mild, except for one of the stars). This section also discusses the derivations of fine- 
structure population ratios, which may be compared to the theoretically expected values 
for different conditions (§6.1). Before one can derive useful limits for the thermal pressures, 
the allowable ranges of temperature must be constrained, and different methods of deriving 
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these constraints are discussed in §6.2. Ultimately, the limits for the population ratios and 
temperatures restrict the possible values for p/k [and local density n(H)], as shown for the 
four cases in the diagrams presented in Fig. 7. Three out of the four stars indicate internal 
thermal pressures for the foreground clouds that are below the generally accepted value for 
the LB. Possible explanations for this imbalance, duplicating that seen for LIC, are presented 
in §7. 

2. Selection of Target Stars 

Target stars chosen for the survey had to satisfy four fundamental criteria. First, the 
stars had to be within the Local Bubble or near its edge. Second, they had to have sufficient 
neutral gas in front so that there was a reasonable expectation of seeing the C 1 absorption 
features. Third, the stars had to be bright enough to give a good signal-to-noise ratio in a 
reasonably short observing time. Finally, the survey avoided stars with projected rotational 
velocities vsini < 50km s~^ (Uesugi & Fukuda 1981 - with actual values provided by 
the VizieR web site at the Strasbourg Data Center), so that stellar features would not 
cause confusion when the interstellar lines were being identified and measured. To satisfy 
the first two requirements, the selection included stars less than about 100 pc away that 
had interstellar D-line absorption features indicating 10^^ < A^(Na 1) < 3 x lO^^cm"^. A 
compilation by Welsh et al (1994) was a good source of information about D-linc absorption 
at the time the survey was being planned. Fluxes at 1565 A listed in the TD-1 Catalogue of 
Stellar Ultraviolet Fluxes (Thompson et al. 1978) provided a good guide for selecting targets 
with satisfactory brightness levels. 

To increase the chance that the C 1 features arose from truly intervening material 
rather than circumstellar disks or shells around the target stars, the survey did not include 
candidates that had spectral types with an emission-hne "e" designation (Hoffleit & Jaschek 
1982; Slettebak 1982) or an excess IRAS flux at 12//m relative to the normal expectation. 
Short-period binaries were also rejected, since they could have interacting gas streams. As a 
flnal precaution against skewing the results with cases dominated by gaseous material very 
near the stars, there was an exclusion of targets that had excess diffuse infrared emission 
nearby in the sky (Gaustad & Van Buren 1993), signifying the possible presence of dust 
grains being heated by the star. 

Four stars that were ultimately chosen for the survey are listed in Tabic 1. Their 
locations with respect to the Local Bubble boundaries mapped by Sfeir et al. (1999) are 
shown in Figure 1. 
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Fig. 1. — Locations of the four program stars listed in Table 1 with respect to the shape 
of the Local Bubble in the plane of the Galaxy, outlined by contours of increasing Na I 
D2-line absorption equivalent widths mapped by Sfeir et al (1999). The Galactic center 
direction is toward the right. To maintain consistency with Sfeir et al's construction of the 
contours according to a spherical projection, the distances of the stars from the origin are 
not foreshortened by cos (6). 
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Table 1. Target Stars and Dataset Identifications 



Star 


HD 


t 


b 




Spectral 


HST Archive 


Name 


Nr. 


(deg) 


(deg) 


(pc) 


Type*" 


Dataset Name 


7 Ori 


35468 


196.9 


-16.0 


75 ±6 


B2III 


Z3CI0107T 


A Lup 


133955 


326.8 


+11.1 


125 + 11 


B3V 


Z3CI0307T 


<^Cyg 


186882 


78.7 


+10.2 


52 + 1 


B9.5IV 


Z3CI0407T 


a Del 


196867 


60.3 


-15.3 


74 + 4 


B9IV 


050001010-40 



'^From the parallax measurements of Hipparcos (Ferryman ct al. 1997) 
reported on the VizieR web site at the Strasbourg Data Center. 

^Hoffleit & Jaschek (1982). 
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3. Observations 

The Goddard High Resolution Spectrograph (GHRS) on HST had a proven capabihty 
of achieving large signal-to-noise ratios for bright targets (Cardelli & Ebbets 1994; Cardelli 
1995). For this reason, this instrument was considered ideal for the objective of recording 
the weak C I absorption features. The survey was originally intended to be completed before 
GHRS was to be replaced by the Space Telescope Imaging Spectrograph (STIS) during the 
second HST servicing mission, but useful observations of a Del were not accomplished until 
after STIS was installed. 

Observations of 7 Ori, A Lup and 5 Cyg were carried out with the Ech-A grating on 
GHRS, with a substep pattern that created 4 spectrum bins per detector diode (STEP- 
PATT=7) plus a samphng of the background level. The standard option for spectrum shifts 
(FP-SPLIT mode) allowed for an elimination of the detector's fixed-pattern noise in the 
spectrum (§4.1 and Appendix A). The GHRS observations were performed with COSTAR 
in place to correct for the telescope's spherical aberration, and the 2" entrance aperture was 
used to increase the efficiency. The STIS observations of a Del with the E140H grating had 
to be done using neutral density filters to prevent the MAMA detector count rate limits 
from being exceeded. The wavelength resolving powers of the GHRS observations were 
A/AA = 90,000 [but with some extended wings in the fine-spread function - see Howk, 
Savage & Fabian (1999)], while that for the STIS spectrum of a Del is 110,000 (Kimble et 
al. 1998; Leitherer 2001). 

The fimited wavelength coverage of the GHRS Digicon detector for echelle spectroscopy 
limited observations to only one multiplet at any particular moment. The multiplet of 
choice was the one at 1260 A, since it had well separated fines and was near the peak in 
the instrument's sensitivity. However, for late B-type stars {5 Cyg and a Del), a strong 
stellar feature seriously depresses the flux level at the position of this multiplet, so other 
multiplcts had to be viewed for these stars. The STIS observation of a Del covered many 
multiplets, but only the strong ones listed in Table 2 gave useful results. For all stars, the 
S II triplet with features at 1250.584, 1253.811 and 1259.519 A was also observed to monitor 
the presence of a generally undepleted element in its favored stage of ionization in an H I 
region. The abundance of S II served a useful purpose in constraining the temperature of 
the gas, as outfined in §6.2. 
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4. Data Reduction 

4.1. Removal of Detector Artifacts 

In order to realize the full potential of GHRS in sensing very weak absorption features, 
it was necessary to remove spurious signal deviations caused by small changes in photo- 
cathode sensitivity with position. The intentional displacements in wavelength for different 
subexposures in the GHRS FP-SPLIT observing routine allowed these variations to be dif- 
ferentiated from real absorption features in the raw spectra. The method used here to derive 
independently the detector's fixed-pattern signal and the real spectral signal gave a direct 
solution and thus differed from the iterative technique described by Ebbets (1992), Cardelli 
& Ebbets (1994) and Fitzpatrick & Spitzer (1994). Details of this more direct method of 
solution are given in Appendix A. 

It is important to recognize that for any given wavelength offset within the FP-SPLlT 
procedure, there are two effects that can alter the mapping of photocathode points onto given 
spectral bins. One of them is a change caused by the Doppler-shift correction for orbital 
motion (performed internally by GHRS). The extremes of this correction can in principle 
differ by as much as twice the orbital velocity of approximately 7 km s~^ times the cosine of 
the angle of the target with respect HST orbital plane at the time of observation. The other 
effect, a much smaller one, arises from small changes in the influence of the Earth's magnetic 
field on the trajectories of the electrons traveling from the detector's photocathode to the 
sensing diodes (Ebbets 1992). Figure 2 illustrates the importance of recognizing these shifts 
and compensating for them in the analysis that eliminates the fixed-pattern noise. The top 
panel shows a straight average of the spectra of 7 Ori (i.e., no attempt to correct for the fixed- 
pattern noise), after shifts had been implemented to compensate for spectral motions with 
respect to the diodes. The middle panel shows the outcome for an initial attempt to correct 
for the fixed-pattern variations by assuming that they remained stationary with respect to 
the pixel assignments in the accumulated signal transmitted to the ground. Practically no 
improvement over the simple average is apparent after this correction. However, one may 
track the movement of the photocathode by measuring the shifts of one or more particularly 
strong fiaws or, if they are not apparent, by cross correlating the spectra. After this is done 
and the spectra have had their fixed-pattern features aligned, an analysis (with appropriate 
modifications in the spectral shifts) produces a satisfactory result (bottom panel). However, 
at the point that we declare the pattern to be fixed with respect to the photocathode, we 
abandon our ability to compensate for differences in response of the detector's diode elements. 
(By an extension of the analysis method presented in Appendix A, it is possible in principle 
to correct for both sources of variation.) 
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Fig. 2. — Different reductions applied to the spectrum of 7 Ori: (top panel) a simple co- 
addition of the individual spectra with shifts that line up the spectral features, (middle panel) 
a composite spectrum corrected for pattern noise assuming that it is fixed with respect to the 
detector diodes, using the method described in Appendix A, and (bottom panel) once again a 
spectrum corrected for the detector pattern, but this time with offsets of up to AA = 0.0348A 
(8.3 km s~^) to compensate for the fact that the strongest ffaws move with respect to the 
diodes from one exposure to the next. See §4.1 for details. The absorption feature on the 
left-hand side is the 1260.736 A transition of C I at a heliocentric velocity of +31 km s~^ 
(uncalibrated for small wavelength errors in the spectrograph). A strong transition of C I* 
at 1261.122 A is too weak to be clearly seen (it should appear at 1261.252 A on the wavelength 
scale presented here). 
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The fixed pattern removal is beneficial only when the fiuctuations are larger than the 
random noise arising from the counting of photoevents. For each spectrum a check was made 
to insure that this condition applied, and, if not, a simple coaddition was performed instead. 
There was only one instance (the coverage of the 1560 A multiplet of C I for 6 Cyg) where 
fluctuations in the fixed-pattern corrected spectrum exceeded those in the simple coadded 
spectrum, and then by only a small amount. 

4.2. Scattered Light Corrections 

Echelle spectrographs create complex background light patterns arising from scattering 
from the echelle and cross-disperser gratings (Cardelli, Ebbets, & Savage 1993). Scattered 
light corrections supplied by the standard calibrated data products supplied by the Space 
Telescope Science Institute were utilized for the GHRS spectra. For the STIS spectrum 
of a Del, the characterization of scattered fight developed by Lindler & Bowers (2000) was 
invoked by utilizing the CALSTIS reduction procedures developed for the STIS Investigation 
Team. Since the C I features toward all of the stars in this survey are weak, errors in the 
scattered fight corrections have a negligible influence on the results. 

4.3. Continuum Definitions 

Least-squares flts of Legendre polynomials to intensities on either side of the absorption 
features deflned the reference continuum levels for the line measurements. The procedures 
of Sembach & Savage (1992) were adopted for determining the appropriate order of the 
polynomial. However, from mock continuum fitting exercises where no interstellar features 
were present, it appeared that terms in the error matrix for the polynomial coefficients 
generally underestimated the true uncertainties in the outcomes by about a factor of two. 
The additional deviations probably arise from errors in the assumptions that polynomials are 
truly appropriate for describing stellar continuum levels over arbitrary wavelength intervals. 
Thus, to account for the uncertainties that are fikely to reach beyond the formal errors, the 
±2(7 deviations were declared to be really ±1(7 uncertainties. 
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5. Analysis of the Absorption Features 

5.1. Column Densities 

Table 2 lists the lines and their measured equivalent widths for the carbon absorption 
features^ in the spectra of each of the four target stars. An example of a case {S Cyg) where 
lines from both C I and C I* could be seen is shown in Fig. 3. 

Errors in equivalent widths arising from the continuum uncertainties discussed above 
(§4.3) were evaluated by comparing the measurement outcomes for the most probable con- 
tinuum levels with those using reasonable departures for the continua (at the declared ilcr 
limits). The resulting uncertainties in Wx should be uncorrelated with errors arising from 
random noise in the signal inside the measurement interval for the line. Hence, in each case 
the estimates for the amplitudes for these two different kinds of errors could be added in 
quadrature. The combined error predictions are reflected in the uncertainties in equivalent 
widths listed in the last column of Table 2. In many cases, the measured values for absorption 
out of C I* were comparable to or weaker than their associated errors, but these numbers 
were retained in the subsequent analysis up to the point where the level populations are 
expressed (in terms of a parameter called /I, as described later in this section), in order to 
preserve the integrity of the final error estimates. 

While the lines are generally very weak, it is a mistake to assume that under all cir- 
cumstances they are either completely unsaturated or fully resolved by the instrument. Very 
high resolution spectra of Na I features observed with ground-based telescopes reveal fea- 
tures with velocity dispersions that are extremely narrow. Since both C I and Na I represent 
ionization stages below those favored for H I regions, their abundances are driven in a very 
similar fashion by local conditions (i.e., primarily the strength of ionizing radiation and local 
electron density). It is thus appropriate to use Na 1 as a surrogate for C 1 when we wish 
to understand how the lines in the present study might saturate. The only complication in 
this comparison arises from the mass difference of the two elements. If turbulence is the 
dominant source of broadening, the profiles of C I and Na I should be virtually identical. 
Conversely, if thermal broadening dominates, the C I hues are broader by a factor (23/12)°-^. 
In the analysis that produced column densities of C I, both extremes were considered. 

Various investigators have determined the nature of the Na I absorptions toward the four 



^Throughout this paper, the designation "C I" appUes to neutral carbon atoms in the ground ^Pq state, 
whereas "C I*" designates the first excited fine-structure level (^Pi). Absorptions from the second excited 
level, C I**=^P2, were too weak to observe. The expression C Itotai refers to neutral carbon in all stages of 
excitation. 
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1277.245^ 


2.225 


2.68 ±0.43 

O.zb ± U.OO 




C I* 


1656.267 


1.987 


1.45 ±0.84 






1657.379 


1.765 


0.43 ± 0.84 






1657.907 


1.890 


1.82 ± 0.73 






1560.682; .709^ 


2.099 


0.73 ± 0.74 






1329.085; .100; . 


123^ 2.166 


0.30 ±0.80 






1277.282^^ 


2.017 


1.11 ±0.45 
0.12 ±0.42 






1277.513*^ 


1.703 


0.52 ±0.56 
-0.14 ±0.50 



^Sources of /-values: multiplets at 1260 and 1329A from Jenkins & Tripp (2001), multiplets at 1560 
and 1657A from Wiese et al. (1996). 

'^Other absorptions from C I* in this multiplet were not measured because they appeared in a small 
wavelength interval where fixed-pattern corrections were not satisfactory enough to reveal weak absorp- 
tion features. See Fig. 2. 

■^Two lines that are unresolved. 

■^There arc two opportunities to measure this feature, since it is visible at opposite ends of two 
adjacent orders of diffraction in the STIS echelle format. 

'^Three lines that are unresolved. The allowance for line saturation was for a feature one-third the 
strength of this ensemble, since the individual features are well separated compared to the real line 
widths and have approximately the same strength. 
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Fig. 3. — Spectrum of S Cyg covering the 1560 A multiplet of C I. The strong absorption 
feature is from the 1560.309 A transition of C I, while the weaker one arises from a blend of 
two features from C I* at 1560.682 and 1560.709 A. 
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stars in this study of C I. Vallerga et al (1993) observed a single absorption component with 
a velocity dispersion b = 0.33 km toward 7 Ori. Crawford (1991) found a best fit to a 
single component in the spectrum of A Lup with b = 1.5 km s~^. Both Blades, Wynne- Jones 

6 Wayte (1980) and Welty, Hobbs & Kulkarni (1994) report very similar values for b for a 
single component in the spectrum of 6 Cyg (the former lists b = 0.44 km s~^ while the latter 
gives b = 0.42 km s~^). The situation for a Del is more complex than the others: Vallerga 
et al (1993) found one component with A^(Na I) = 9.5 x 10^° cm-^ and b = 1.10 km s~^ and 
another, slightly overlapping feature with A^(Na I) = 1.6 x lO^'^cm"^ and b — 0.90 km s""^. 

Table 3 lists the column densities of C 1 and C 1* based on curves of growth for the 
two extreme assumptions for the character of the line broadening: (1) the line shapes are 
exclusively determined by thermal Doppler motions at a temperature 

_ 23mp6^3^ J 

2k ^ ' 

so that 6c I = (23/12)°'^6Na ii and (2) the temperatures are arbitrarily small and the broad- 
ening arises principally from some form of bulk motion (e.g., turbulence), making the 6- values 
for the two elements identical. For either of the two cases, the curvc-of-growth corrections 
were explicitly calculated for each line using the nominal W\ values and their accompanying 
±1(7 error limits. When more than one transition could be used (or a single transition was 
viewed more than once - see Table 2) , the results reflect a weighted average, with weights 
proportional to cr"^ and a final error equal to (X^j c~^)~'^'^. For the strongest C I lines, 
typical saturation corrections arising from the curves of growth were as follows: 0.05 dex for 

7 Ori, 0.01 dex for A Lup, 0.27 dex (thermal) or 0.51 dex (turbulent) for 8 Cyg, and 0.08 dex 
for a Del. 



5.2. Fine-Structure Excitations 

A conventional means of expressing the fine-structure excitation of C I is through the 

quantities /I = N{C r)/N{C Itotai) and /2 = A^(C r*)/N{C Itotai) (Jenkins & Shaya 1979; 
Jenkins et al. 1981, 1998; Jenkins, Jura, & Loewenstein 1983; Smith et al. 1991; Jenkins 
& Wallerstein 1995). The utility of this representation arises from the ease of interpreting 
discrepancies in the excitations of the two levels in terms of a superposition of contributions 
from regions with contrasting physical conditions (Jenkins & Tripp 2001). Since C 1* was the 
only excited level that created strong enough lines to measure the present survey, the special 
advantage of the /l,/2 representation is lost. Nevertheless, in the interest of maintaining 
consistency with the earlier studies of conditions elsewhere, there is some merit in retaining 
the /I representation in favor of the more straightforward quantity A^(C I*)/A^(C I). 



-14- 



Table 3. Column Densities and /I Values 



Star Line Assumed T N{C I) N{C I*) fl 

Name Broadening'' (K) (lO^^cm"^) (lO^^cm"^) 



7 Ori 


thermal 


152 


1.13 ±0.12 


0.104 ±0.171 






turbulent 


< 152 


1.17±0.14 


0.105 ±0.171 




A Lup 


thermal 


5560 


1.10±0.16 


0.251 ±0.177 


0.1791°:}^ 




turbulent 


< 5560 


1.11±0.17 


0.251 ±0.177 


n 1 70+0-122 
'^■^ ' ^-0.125 


SCyg 


thermal 


245 




1.27 ±0.40 


U.iO4_Q057 




turbulent 


< 245 




1.36 ±0.45 


A -,Ao+0.063 


a Del 


thermal 


1680 


1.91 ±0.13 


0.51 ±0.18 


A 204+0-064 




turbulent 


< 1680 


2.07 ±0.14 


0.51 ±0.18 


A -| 00+0.061 
'-'■-'-^■^-0.065 



^Calculations were performed for two possible extremes in the assumed line 
broadening of the C I features: either the lines broadened by thermal Doppler 
motions for a temperature given by Eq. 1 (sec the numerical values in the next 
column) or they are broadened purely by turbulence with arbitrarily low gas 
temperatures. 

^/l = A^(C I*)/A^(C Itotai)- Absorption features from the second excited 
level, A^(C I**), arc not observable because they are too weak. Hence, a value 
of N{C I**)/iV(C I*)=0.15 is assumed to arrive at a value for N{C Itotai)- See 
text for details. 
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To evaluate /I, one needs to estimate the amount of C I** present, since it makes a 
small contribution to C Itotai- Calculations of the expected level populations indicate that 
A^(C I**) ~ 0.15A^(C I*) when the excitations are low and there is no complex mixture of 
gases with very different conditions [e.g., see the theoretical tracks in Fig. 6 of Jenkins & 
Tripp (2001)]. Thus, values of fl listed in the last column of Table 3 were evaluated from 
the expression 

n ^ (2) 

■' N{C I) + 1.15A^(C I*) ■ ^ ^ 

Upper or lower limits for fl represent the outcomes from Eq. 2 using opposite 1-a extremes 
of the permitted column densities of A'"(C I) and A^(C I*). 



6. Interpretation 

6.1. Expected Values of fl 

The next step in the investigation is to explore the expected behavior of fl over a broad 
range of possible physical conditions. The equilibrium concentrations of the atoms in dif- 
ferent fine-structure excitations are governed by the balance between collisional excitations, 
coUisional de-excitations, and the levels' rates of spontaneous radiative decay. Small adjust- 
ments arise also from optical pumping of the levels (Jenkins & Shaya 1979). The calculations 
adopted here duplicate those described by Jenkins & Tripp (2001), except for an additional 
task of estimating the fractional ionization of the gas in the regime of high temperatures 
and low densities. Knowing the degree of ionization is important because electrons and pro- 
tons have collision rate constants that are appreciably higher than those of neutral hydrogen 
atoms when temperatures are much above 100 K (Keenan 1989). 

The diffuse, neutral interstellar medium is partly ionized by cosmic rays (at an estimated 
rate ^cr = 3 x lO^^^s"^ that includes ionizations from secondary particles) and EUV and 
soft x-ray photons (at a rate ^xr) that manage to penetrate through an absorbing layer of 
neutral H and He near the edges of a region (Wolfire et al. 1995). The fractional ionization 
of H is governed by the balance between these ionizing processes and recombinations, 

n(H I)(^cR + ^xr) = Q;2,Hri(e)^ + ag^iin{H)n{e) (3) 

where the right-hand side of the equation includes two recombination channels for the pro- 
tons: (1) recombinations with free electrons to all energy states n — 2 and higher and (2) 
recombinations with negatively charged, small grains and polycyclic aromatic hydrocarbon 
molecules (Weingartner & Draine 2001), with a rate constant Ug^n normahzed to the total hy- 
drogen density n{H). If we assume the shielding of the inner region of a cloud arises from an 
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external column density A'"(H I) ~ 10^^ cm~^ (an approximately correct value for the regions 
viewed here; see §6.2), the calculations of Wolfire et al (1995) indicate that ^xr = 7x 10~^^ 
for n(e) ~ 10~^ cm~^. Supplementing the free electrons from the ionization of H (and He) are 
those created from elements that are almost completely ionized because they have ionization 
potentials less than 13.5 eV. This contribution is about 2 x 10~*^n(H). 

Figure 4 shows the outcomes for the expected values of /I in the form of a contour 
diagram in a representation of \og{p/k) vs. logn(H). It is clear that without independent 
information that constrains either T or n(H), measurements of /I are of little value in defining 
the thermal pressure p (except for defining lower bounds when there is a well defined lower 
limit for /I). In the following discussions, we consider different ways to constrain the allowed 
temperatures in the C I-bearing regions. Even rudimentary limits can be beneficial. 

6.2. Temperature Ranges 

6.2.1. Upper Limits 

The most straightforward upper limits for T arise from the narrowness of the Na I lines 
seen in the spectra of 7 Ori and 5 Cyg. The application of Eq. 1 to the 6-values given in 
§5.1 yields upper limits T — 152 and 245 K for the respective stars. These temperatures are 
located only slightly above the sequence of positions where the fl contours loop down to 
their minima in \og{p/k) in Fig. 4. Unfortunately, the broader lines for Na 1 toward A Lup 
and a Del could allow much higher values of T, i.e., ones that do not sufficiently constrain 
the allowed p/k much below 10^'^cm~^K for the measured values of fl toward these stars. 
Additional means for constraining the allowed temperatures must rely on some independent 
theoretical arguments. 

The tactic for determining the largest probable temperatures of material in front of 
A Lup and a Del will be to prove that above a certain temperature, it is difficult to produce 
the observed amount of C I. In order to do this, we must start with an estimate of the 
amount of gas present in each of the C I-bearing clouds. It is generally acknowledged that 
the abundance of sulfur atoms in the interstellar medium is close to the solar abundance 
ratio relative to hydrogen, S/H = 1.9 x 10~^ (Federman et al. 1993; Spitzer & Fitzpatrick 
1993; Fitzpatrick & Spitzer 1994, 1997; Howk, Savage, & Fabian 1999), i.e., sulfur is not 
appreciably depleted onto dust grains. Sulfur in its singly ionized form is the most abundant 
stage expected for H 1 regions. Since S II can also reside in H II regions (its ionization 
potential equals 23.4 eV) a determination of A^(H I) using A^(S II) is, strictly speaking, 
only an upper limit. However the arguments about not producing enough C Itotai to match 
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Fig. 4. — Expected values of /I = A^(C I*)/A^(C Itotai) [solid contour levels] as a function 
of the parameters \og{p/k) and logn(H). Diagonal straight, dotted lines show the loci of 
constant temperatures T with labels indicating the values. Curved, dashed lines show the 
thermal equilibrium curves derived by Wolfire et al. (1995) for two different amounts of 
shielding of EUV and soft x-ray radiation (upper curve: A^(H I) = 10^^ cm^^; lower curve: 
7V(H I) = 10^9 cm-2). 
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the observations will only be strengthened if A^(H I) is really less than this limit. Thus we 
can simplify the discussions and adopt a conservative position by treating the estimate for 
A^(H I) as an actual value. 

The middle and lowest panels of Figure 5 show the two weakest absorption features of 
S II in the spectra of a. Del and A Lup. The ratios of the /-values of the two transitions is 2.0 
(Morton 1991), but it is clear from the plots that the ratios of the apparent optical depths in 
the stronger parts of the profiles are smaller than 2. While one way to determine A^(S II) is 
to measure the equivalent widths of the lines and use the standard curve of growth analysis, 
a better approach is to measure the apparent optical depth ratios of the two lines at every 
velocity and then apply a method described by Jenkins (1996) to correct for distortions in 
the weak line caused by unresolved, saturated structures within the profile. 

For a Del, A^(S II) = 3.1 x lO^^cm-^. In turn, one infers that iV(H I) = 1.6 x lO^^cm-^ 
and 7V(C II) = 2.2 x 10^^ cm~^ if the abundance of carbon relative to hydrogen is the typical 
value of about 1.4 x 10"'^ in the interstellar medium (Hobbs, York, h Oegerle 1982; CardeUi 
et al. 1991, 1993, 1996; Sofia et al. 1997; Sofia, Fitzpatrick, & Meyer 1998). Since the 
ionization and recombination properties of S II and C II are very similar (Sofia & Jenkins 
1998; Howk & Sembach 1999), shifts in the ratio of these two first ions caused by ionization 
in low density, partly ionized gas should not be a problem. Comparing the neutral carbon 
column densities toward a Del reported in Table 3 to the inferred A^(C II) derived here, 
we arrive at logA^(C Itotai) ^ logA^(C II) = —2.97. The significance of this number will be 
evident when it is compared to theoretical expectations under many different conditions. 

Figure 6 shows the contours for the expected values of log n(C Itotai) — log n(C II) arising 
from the equilibrium condition 

n(C ItotaO^c I = n(C II) [n(e)ac + a3,cn(H)] (4) 

where the assumed interstellar ionization rate i = 2.24 x 10^^*^ s^^ (Jenkins & Shaya 1979), 
the recombination rate of singly-ionized carbon with free electrons etc is evaluated from the 
fitting equation given by ShuU & van Steenberg (1982), and the grain/PAH recombination 
rate is from the relation specified by Weingartner & Draine (2001) with their radiation 
field strength parameter G set to 1.13. Electron densities n(e) arise from solutions to Eq. 3, 
as described earlier. [If one neglects grain recombination altogether, the predicted values for 
n(C Itotai)/'T'(C II) are not much different because the reduced total recombination rate is 
compensated by a larger predicted value for n{e)\ 

By comparing the contours in Figs. 4 and 6, it is clear that if a. Del were to have /I 
equal to the nominal value of 0.20 (for pure thermal broadening of the lines), the temperature 
must be equal to or less than 150 K to give an expected yield of C Itotai that is at least as 
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Fig. 5. — Absorption profiles for various transitions of C I and S II, as indicated, in the 
spectra of a Del (left) and A Lup (right) . The intensities have been normalized to the stellar 
continua. Note that the scale for the C I transition for A Lup has been expanded. 
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Fig. 6. — Expected values of logn(C Itotai) — logn(C II) as a function of the parameters 
log(p/A;) and logn(H), found by solving Eq. 4. Other lines in the diagram are described in 
the caption for Fig. 4. 
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large as the observed value. If /I is as high as the +la limit, the temperature could be as 
high as 250 K. The upper limit for T if fl is at its —la limit is 90 K. These temperature 
limits are substantially lower than those set by the width of the Na I absorption feature. 

For A Lup, the strong, central peak of the S II absorption matches the velocity of the 
C I profile (see Fig. 5). The S II associated with this peak has a column density iV(S II) = 
1.2 X 10^^ cm^^. (It is reasonable to consider the gas giving absorptions at velocities on either 
side of this peak to arise from unrelated material.) Using reasoning identical to that applied 
to a Del, one finds for A Lup that iV(H 1) = 6.2 x lO^^ cm'^, A^(C II) = 8.7 x lO^"' cm'^, 
and log A^(C Itotai) ~ log A^(C II) = —3.81. Clearly, the restriction on T for A Lup is not as 
stringent as for a Del: a limit T < 2700 K arises from the interception of a line corresponding 
to fl — 0.19 (nominal value) with the expectation logiV(C Itotai) — logA^(C II) = —3.81. 
This limit is only mildly more restrictive than the one that arises from the line width of 
Na I. For fl — 0.05 (the —la hmit) the limiting temperature drops to 350 K. 

6.2.2. Lower Limits 

A task which now remains is to confine the allowable temperatures on the low side, since, 
as with high temperatures, arbitrarily high pressures could occur if T were permitted to be 
arbitrarily low (again, see Fig. 4). In order to realize such a constraint, one may consider 
the sequence of temperatures at different densities predicted for thermal equilibrium in the 
interstellar medium. Wolfire et al. (1995) calculated equilibria for the diffuse phases of 
interstellar material for various assumptions. The two curved, dashed lines in Fig. 4 show 
the outcomes for their predictions based on heating rates for material shielded by absorbing 
columns A^(II I) = 10^^ (upper curve) and lO^^cm"^ (lower curve), using cooling rates 
based on normal grain and heavy element abundances. According to Field's (1965) thermal 
stabihty criterion, stable phases of the interstellar medium should be located only along 
portions of the equilibrium curves that have a positive slope in the logp vs. logn diagram 
(Field, Goldsmith, & Habing 1969; ShuU 1987; Begelman 1990). In a simplistic apphcation 
of this principle, one would predict that stable phases would exist only at shghtly below 
10,000 K (left-hand branch) or in the range 30 < T < 300 K (right-hand branch). In reality, 
this picture is only partly correct. On the one hand, observations of 21-cm emission and 
absorption in directions toward extragalactic sources indicate that about half of the H I in the 
diffuse medium violates this condition (Heiles 2001). The data show strong evidence that gas 
can often be found at temperatures that are intermediate between the two (positive-slope) 
equilibrium curves, an effect that is probably attributable to the effects of mixing of the two 
extremes in a regime strongly dominated by turbulence (Vazquez-Semadeni, Gazol, & Scalo 
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2000; Gazol et al. 2001). On the other hand, the observations generally show that the lowest 
H I temperatures coincide with the theoretical expectations. While broad surveys of emission 
with superposed absorption highlight a few, selected locations where spin temperatures are 
below 30 K (Gibson et al. 2000; Knee & Brunt 2001), the sight lines toward extragalactic 
sources offer a chance to sample random volumes: the survey of emission and absorption 
toward these sources reported by Heiles (2001) and Heiles & Troland (2002) reveals that only 
about 3.7% of the cold, neutral medium is at temperatures below 25 K. In short, the 21-cm 
data indicate that a good working assumption is that T is usually above the right-hand 
branch of the thermal equilibrium curve, and violations of the lower limit defined by this 
curve are rare. For this reason, it is reasonable to regard the nearly coincident dashed lines 
on the right-hand side of Fig. 4 as a good representation for the lowest probable temperatures 
for different pressures. 



7. Conclusions and Possible Interpretations 

A synthesis of the conclusions presented earlier is shown in Figure 7, again in the 
representation of log(p/fc) vs. logn(H). To constrain the allowed physical conditions, the 
outcomes presented in Table 3 for the fine-structure population ratios of C 1, represented 
by the parameter /I, must be supplemented by intersecting lines representing temperature 
limits derived in §6.2. Of these limits, the upper bounds derived from the measurements of 
the Na 1 line widths by other investigators are the most direct and reliable. However these 
measurements are valuable only for the cases of 7 Ori and 5 Cyg. Alternative means for 
limiting T, such as using the thermal equilibrium curve to define a lower bound and the 
carbon ionization equilibrium for an upper bound, are indirect and subject to less certain 
theoretical assumptions. 

It should be noted that for each case the lines that border the allowed (shaded) regions 
do not enclose the worst possible deviations. For instance, the /I limits represent only ±1(7 
deviations. The temperature constraints from the width of the Na I absorption lines do not 
include possible errors in measurement, since they are difficult to assess. The theoretical 
limits based on thermal equilibria or the expected fractional abundances of carbon in neu- 
tral form are only as good as the assumptions that were incorporated into their respective 
developments. In spite of these shortcomings, when three of the stars (those other than 
A Lup) are examined collectively, they indicate that clouds within the Local Bubble have 
thermal pressures in the range 10^ < p/k < 10'^ cm K. 

From the introductory remarks presented in §1, it is clear that the most interesting 
limit on pressure is the upper one. Any single case is not compelling enough to support the 
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log n(H) 

Fig. 7. — Allowed combinations of log(p/A;) and log[n(H)] permitted by the ±1(7 error limits 
for /I for the C I-absorbing regions in front of the four targets. Intersecting constraints arise 
from the thermal equilibrium curve for an external absorbing column A^(H 1) = 10^^ cm~^ 
calculated by Wolfirc, ct al (1995), and temperature limits defined either from the Doppler 
width of Na 1 absorption lines or, alternatively, the requirement that the expected yield of 
C Itotai relative to C 11 from the condition of ionization equilibrium (Eq. 4) is as least as 
much as actually observed. General layouts of the expected C I fine-structure population 
ratios and C Itotai/C II ionization equilibria are shown in Figs. 4 and 6, respectively. 
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proposal that p/k < lO^cm"^ K, since deviations in excess of Icr in a particular direction can 
be expected to occur 16% of the time. Nevertheless, if one supposes that all clouds within 
the Local Bubble have a single representative pressure, it is difficult to imagine that all of the 
random observational errors would conspire to mislead us to believe that p/k < 10^cm~^K 
when in reality the pressure is much higher than this value. The result for A Lup is not 
helpful in constraining p on the high side, but at least it shows an outcome on the low side 
that is consistent with those of the other stars. As for systematic errors, the only easily 
identifiable ones are possible errors in the relative /-values of the C I and CI* transitions or 
errors in the parameters used to develop the predicted values of /I, such as the coUisional 
rate constants or spontaneous radiative decay rates. 

It is apparent from this study that a seemingly large disparity in pressure between 
neutral clouds embedded in the Local Bubble and the surrounding hot gas is a general phe- 
nomenon and not just confined to the LIC around the Sun. Possible solutions to the dilemma 
posed by the apparent pressure imbalance may arise from the following considerations: 

1. Perhaps the clouds have not had enough time to respond fully to a sudden increase 
in external pressure that happened during the event (s) that created the LB, perhaps 
some 10^ — 10^ yr ago (Cox & Reynolds 1987). If this is the case, each cloud is in 
an intermediate condition where a shock is still progressing from the periphery of the 
cloud, perhaps well outside the portion containing the C I, to its center. However, 
this explanation is not satisfactory for the LIC, since the distance from the Sun to the 
closest boundary divided by the sound speed is less than 10^ yr. 

2. Perhaps magnetic fields within the clouds provide the support needed to counteract 
the high pressure from the outside. An upper limit B < 3/xG for the magnetic field just 
outside the hehosphere (Gloeckler, Fisk, & Geiss 1997) diminishes the attractiveness 
of this explanation in the case of the LIC. 

3. Additional contributions to the clouds' internal pressures might arise from turbulent 
motions. However, the narrowness of the Na I lines observed for 7 Ori and 5 Cyg lowers 
the prospects for this proposal. For instance, if one uses the limit n(H) < lO^'^cm"^ 
shown in Fig. 7 and assumes the Na I lines are broadened only by turbulence, one finds 
that Pturb/^ — ph^ I'^k < 470 and 800cm~"^K for 7 Ori and b Cyg, respectively. 

4. Perhaps the measurements of the thermal pressure for the hot gas in the Local Bubble 
were incorrect because they relied on either inaccurate models for the fine emission or 
the assumption that the radiation from the hot plasma is emitted under equilibrium 
conditions. For a plasma that is cooling very rapidly, perhaps through adiabatic ex- 
pansion, the x-rays may come principally from delayed recombinations in cooler gas. 
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rather than coUisional excitations at higher temperatures (Breitschwerdt & Schmutzler 
1994). Apart from this theoretical proposal, there have also been persistent difficul- 
ties in reconciling measurements of x-ray and EUV line emission with models for a 
hot plasma in collisional equilibrium and with cosmic abundance ratios of the heavier 
elements (Jelinsky, Vallerga, & Edelstein 1995; McCammon et al. 2002). In trying 
to understand the soft x-ray spectrum measured by the Diffuse X-ray Spectrometer 
(DXS), Sanders et al. (2001) even explored models with depleted element abundances 
and nonequilibrium conditions arising from either rapid heating or cooling. Prom their 
inability to achieve reasonable fits of these more complex models with their data, they 
suggested that the discrepancies might arise from inaccurate atomic data which were 
used to predict the expected emission line strengths. If these models are in need of re- 
vision, then one might question the earlier estimates of hot gas emission measures from 
low resolution surveys of x-ray emission that led to the finding that p/k > 10'^cm~^K 
within the LB. 



The author thanks J. Weingartner for supplying a program that calculated the partial 
ionization of the ISM, as described in §6.1. Important parameters for the target stars were 
obtained from the VizieR web site (Ochsenbein, Bauer, & Marcout 2000), as noted in the text 
and Table 1. Support for this research was provided by NASA through grant number GO- 
6415.01-95A from the Space Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 

A. Fixed-Pattern Noise Analysis 
A.l. Equations 

The discussion here addresses the problem of disentangling the detector's fixed pattern 
response function^ from the spectrum that moves, on the basis of information presented 
in the separate exposures. The analysis method solves for the two functions, but with 



^As noted in §4.1, for GHRS the fine-scale sensitivity pattern is not always stationary with respect to 
the readout diode array. One must explicitly determine the small shifts in position of the photocathode's 
pattern function for each exposure, since the mapping of the photocathode onto the readout diodes can 
change with time. In the discussion presented here, the term "stationary" means stationary with respect 
to the photocathode's reference frame, which may or may not move across the diodes from one exposure to 
the next, depending on the observing circumstances. The notion of a movement refers to the intentional 
displacement of the spectrum caused by a commanded motion of the grating carousel. 
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the actual work being carried out in the Fourier Transform domain. This method has two 
important features. First, it offers a direct evahiation of the solutions, rather than relying on 
an initial guess for the two functions that is followed by iterations that gradually improve the 
representations until they converge toward consistent representations. Second, by exploring 
the properties of equations that operate in the Fourier domain, we gain insights on some near 
singular conditions that, in some practical circumstances, give poorly defined results when 
the shift dimensions have common factors in wavelength. This problem, discussed in §A.2, 
warns us of some potential dangers of programming shifts that are of uniform magnitude, as 
was intended for the FP-SPLIT routines on GHRS. 

An individual spectrum F[X) recorded by a spectrograph consists of the actual spec- 
trum A{X) multiplied by the pattern function P(A), plus a noise contribution arising from 
statistical fluctuations in the recorded photoevents. For the Digicon detector on GHRS, 
P(A) has the appearance of a function that, over large scales, is nearly constant, but with 
a high frequency graniilarity having an amplitude of a few percent. However, occasionally 
there are flaws in the photocathodc that give much larger spikes (Ebbets 1992). In the refer- 
ence frame of the detector's photocathodc, separate recordings of Fi{X) using the FP-SPLlT 
routine consist of the stationary pattern function P(A) with different superpositions of the 
actual spectrum A{X) having offsets with magnitudes 5i in the wavelength direction, where 
the subscript i differentiates the separate exposures. The magnitudes of Si may be deter- 
mined by either cross correlating the exposures (for a spectrum with many narrow features) 
or comparing the offsets of positions of a strong line. The latter is recommended over the 
former if a single, narrow feature is present, with the rest of the spectrum being mostly 
random noise or much broader features. 

As a first step in the analysis, we can reformulate the observed spectrum F{X) in terms 
of its logarithm, so that the mixing of two functions represented by -P(A) and A{X) is a 
simple hnear sum, log[P(A)] -|- log[74(A — Si)]. A restriction that should be imposed is that 
downward excursions of both P{X) and A{X) are not very large compared with the average 
signal levels. If this condition is violated, the negative dips can have a disproportionately 
large influence on the outcome, especially because the analysis works in the Fourier domain 
where large (negative) 5-functions can have a global effect. For P(A) this restriction to small 
amphtudes is generally satisfied for GHRS, except at the locations of serious fiaws in the 
photodiode array (which are normally corrected out beforehand). For the actual spectrum 
A{X), we require that there are no very deep absorption features. (Methods for overcoming 
problems with deep features in the spectrum will be discussed very briefiy in §A.2.) 

In the Fourier domain we define the real and imaginary parts of the pattern's transform 
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by the functions in inverse wavelength a;, 

p,H = 3?{^[logP(A)]} (Al) 

and 

p,M=9{^[logP(A)]}, (A2) 

where T is the Fourier transform operator. For F(A) and ^(A) we designate the functions 
fxi fy, Ojx and Qy in the same manner and allow for the fact that separate recordings will have 
different phase shifts for a that are invoked by multiplying the zero-shift, complex function 
ax{uj) + iay{uj) by a factor exp(— 27ria;(5i). 

The best solutions for Pxioj), Pyi^^), dxi^), and ay{u!) may be evaluated at each cu 
independently by solving for the smallest possible sum of the squared real and imaginary 
residuals 

n 

Ql = ^i(^x cos Ai + ay sin Ai+p^- fx,i)'^ (A3) 

i=l 

and 

n 

Ql = ^(-a^ sin Ai + Uy cos Ai+py - fy^if , (A4) 

i=l 

where Aj is an abbreviation for 27ra;5j. This minimum is achieved by setting the partial 
differentials of with respect to the four unknowns Ux, ay, Px, and Py equal to zero. If 
u(a;) represents a vector of these respective unknowns at each cu, we must solve a system of 
4 linear equations 

C(a;) • u(a;) = b(a;) (A5) 
where the coefficients of Cij of the matrix C and bi of the vector b are hsted in Table 4. 

After solving for each vector u{uj) over all u, one simply evaluates the inverse Fourier 
transforms of the two pairs of terms, Ui{lj) + iu2{uj) and Uz{ijj) + iUi^{ijj) to recover the best 
representations of A{\) and -P(A), respectively. 

A. 2. Limitations 

Large increases in the magnitudes of terms in the error matrices [C(c(;)]~^ at certain 
frequencies indicate places where the solutions are ill-defined. As one would expect, the 
outcomes for frequencies whose inverses are much larger than the total span of the 5, are 
not well determined. The solutions for A{\) and -P(A) can exhibit very broad, spurious 
undulations (of opposite sign), as an outcome of the analysis program's attempt to reconcile 
small differences in the low-frequency components in the observed spectra arising simply 
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from noise and subtle systematic errors. As a practical matter, it is reasonable to assign 
artificially a condition that any low-frequency components of the observed spectra must arise 
from either the source spectrum A{X) or the detector response -P(A). A similar quandary 
arises for certain frequencies where integral numbers of sinusoidal variations in the spectrum 
can exactly match the spacings between all of the 6i. One consequence of this phenomenon 
is that the programmed intent of having equal step sizes the FP-SPLIT routine for GHRS 
was a misguided one.^ Fortunately, in practice the step sizes usually turn out to be not quite 
equal, and this removes the degeneracies of the solutions. 

The presence of deep absorption features in a spectrum poses special difficulties. At 
the locations of these features, the logarithms of F{X) surge to very large negative values, 
becoming even uncomputable in the bottoms of lines where the absorption is total. One 
must acknowledge that there is no way to determine P{X) at locations where less than two 
of the spectra have non-zero fluxes. Inelegant but nevertheless workable methods of avoiding 
irregularities in the computations include adding a small offset to the spectrum so that zero 
intensities are never encountered or, alternatively, eradicating the line and replacing it with 
a continuum intensity level. 



^This shortcoming was noted by the STIS Instrument Definition Team. As a result, the separation of the 
multiple slits for the FP-SPLIT option in STIS are incommensurate with each other (Leitherer 2001). 
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Table 4. Coefficients^ for the Matrix and Vector in Equation A5 



Term 










i 




i = 2 J = 3 







dax 
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n 





E cos Aj 


- X sin Aj 


f^^i cos Ai - fy,' 


i sin Aj 


dgi 

day 


2 
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E sin Aj 


COS Aj 


E/x,iSinAi + E/y,^ 


; cos Aj 


dgi 
dpx 


3 


cos Aj 


sin Aj 


n 









dpy 


4 


- ^ sin Aj 


S cos Aj 





n 


Y fy,i 





'^Each sum is evaluated from i = 1 to n, where n is the number of independent spectral 
recordings. 
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